Differential Correction Map for GNSS

ABSTRACT

A method comprises generating for a Global Navigation Satellite System (GNSS) a differential correction map (DCM) representing a non-planar surface of differential corrections that varies across a geographical area represented by the DCM, the differential corrections being based on a reference station constellation of GNSS reference stations having respective locations spanning the geographical area, the GNSS reference stations including physical or virtual reference stations.

PRIORITY CLAIM

This application claims priority to U.S. Provisional Application No.62/688,260, filed Jun. 21, 2018, which is incorporated herein byreference in its entirety.

TECHNICAL FIELD

The present disclosure relates to differential corrections for GlobalNavigation Satellite Systems (GNSSs).

BACKGROUND

In a conventional networked Differential Global Positioning System (GPS)(DGPS) service, a mobile platform or rover determines its GPS locationand reports the GPS location to the DGPS service. The DGPS servicecomputes a differential correction for a single virtual referencestation at the GPS location, based on information from actual physicalreference stations. The DGPS service transmits the differentialcorrection based on the single virtual reference station back to therover, which uses the differential correction to correct the GPSlocation. The networked DGPS service disadvantageously requirescontinuous two-way communications with the rover. Such two-waycommunications may not be available at all times and locations (coveragearea), and a user of the rover is likely to incur a subscription cost(either directly or indirectly). There is also a potential privacy issuesince the DGPS service learns the location of the rover each time therover reports its location.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is an illustration of an example environment in which local DGPSreference (DGPSR) stations correct GPS errors to achieve precision GPSmeasurements at local vehicles or rovers.

FIG. 2 is an illustration of an example environment for networked DGPSapplications covering a wide area.

FIG. 3A is an illustration of an example environment that supports abroadcast method of delivering DGPS corrections to rovers via adifferential correction map (DCM).

FIG. 3B is an illustration of an example environment that supports abroadcast method of delivering DGPS corrections to rovers acrossmultiple broadcast areas via multiple DCMs.

FIG. 3C is an illustration of an example environment in which DCMs aremade accessible over a communication network.

FIG. 4 is an illustration of example virtual DGPSR stations that form aVirtual Reference Station Constellation (VRSC).

FIG. 5A is an example isometric DCM plot for interpolated DGPSRdifferential corrections for the VRSC of FIG. 4, for exponent p=1 of aproximity-weighted interpolation based on separation distance.

FIG. 5B is an example isometric DCM plot for interpolated DGPSRdifferential corrections for the VRSC of FIG. 4, for exponent p=² of aproximity-weighted interpolation based on separation distance.

FIG. 6A is an illustration of example complementary VRSCs including theVRSC of FIG. 4.

FIG. 6B is an example isometric DCM plot of proximity-weightedinterpolated differential corrections averaged across the complementaryVRSCs of FIG. 6A.

FIG. 7A is an example isometric DCM plot from a quadratic polynomial fitof virtual DGPSR stations.

FIG. 7B is an example isometric DCM plot from a modified quadraticpolynomial fit of the virtual DGPSR stations from FIG. 7A.

FIG. 8 is an example isometric DCM plot from an average fit of modifiedpolynomial and proximity-weighted interpolation methods.

FIG. 9 is an illustration of an example DCM Type 0 forproximity-weighted interpolation based on the VRSC of FIG. 4.

FIG. 10 is an illustration of an example DCM Type 1 based on apolynomial fit for the VRSC of FIG. 4.

FIG. 11 is a flowchart of an example method of generating anddisseminating a DCM that may be performed in any of the environments ofFIGS. 3A-3C.

FIG. 12 is a flowchart of an example method of generating andbroadcasting multiple DCMs across multiple broadcast coverage areasperformed in the network environment of FIG. 3B.

FIG. 13A is a flowchart of an example method of processing a DCM,performed by a processor of a rover.

FIG. 13B is a flowchart of an example method of processing multiple DCMsbased on the method of FIG. 13A.

FIG. 13C is a flowchart of another example method of processing multipleDCMs based on the method of FIG. 13A.

FIG. 13D is a flowchart of another example method of processing a DCMbased on the method of FIG. 13A.

FIG. 14 is a block diagram of an example DCM DGPS management entity usedto generate and assist with broadcasting a DCM.

FIG. 15 is an example block diagram of a rover to receive and processthe DCM.

FIG. 16 is a block diagram of an example radio station to broadcast theDCM.

DESCRIPTION OF EXAMPLE EMBODIMENTS Example Embodiments INTRODUCTION

Global Navigation Satellite Systems (GNSSs) have become crucialcomponents of the quickly-evolving autonomous vehicle navigationtechnology. The location accuracy of ubiquitous consumer-grade GlobalPositioning System (GPS) receivers is typically within about 3 meterswhen 4 or more GPS satellites are in line of sight (LOS). This accuracyis generally sufficient for most vehicle navigation where a GPS-enableddevice is used to provide directions to assist a human driver. In thiscase, the human driver is able to perform the vehicle operations thatrequire finer location precision, such as staying in the lane orparking. However, autonomous vehicle navigation requires greaterlocation precision provided by a combination of devices including GPS.Redundancy is needed since the vehicle sensors are not each reliable allof the time. For example, the vehicle's GPS antenna may be shadowed fromthe satellites by buildings or terrain, and multipath signal reflectionscould cause errors in GPS measurements. Improvements in accuracy beyondthat of typical consumer-grade GPS devices require a variety oftechniques to yield GPS location accuracy within about 10 cm in LOSconditions, sufficient for autonomous vehicle navigation.

A GPS receiver measures the relative time delays of the signals frommultiple GPS satellites to compute its location, time, and velocity. Theerrors in these time-delay measurements are due to multiple factors. TheGPS signals from the satellites are subject to satellite clock andephemeris (orbit trajectory) errors. The LOS signals arriving at the GPSreceiver antenna experience errors due to ionosphere and tropospheresignal propagation (index of refraction variations). The GPS receiverincurs measurement errors due to imprecise signal processing algorithmsas well as unwanted signal reflections from objects. Redundant devicessuch as Inertial Navigation Systems (INS), image processing, radar andsonar are used to aid local navigation and estimate location when theGPS signals are not available or faulty (multipath reflections).

The GPS measurement accuracy over single frequency (L1) consumer-gradeGPS devices can be improved by a combination of techniques. Dualfrequency (L1 and L2) GPS receivers offer greater inherent measurementaccuracy due to the availability of more GPS satellite signals.Real-Time-Kinematic (RTK) techniques use not only the demodulated GPSpseudo-noise (PN) signal for relative time measurements, but also usethe carrier frequency cycles for much finer resolution. These receiverprocessing improvements alone would not yield their inherent accuracypotential unless the other errors of the GPS signal are also corrected.This is accomplished with a class of techniques called Differential GPS(DGPS), where a precision GPS receiver accepts a DGPS correction input(e.g., in Radio Technical Commission for Maritime Services (RTCM)format).

DGPS Reference Station

DGPS techniques provide auxiliary error correction information to theGPS receiver. This error (differential) correction information about thelocal GPS signals includes ephemeris and clock errors from eachsatellite in view, as well as error correction for local ionosphere andtroposphere propagation. The simplest method of providing DGPScorrection data is accomplished by measuring the differential error at alocal DGPS reference (DGPSR) station. The DGPSR station is an actual,physical DGPSR station, also referred to as an installed DGPSR station.

FIG. 1 shows an environment in which local DGPSR stations 102 (referredto collectively as DGPSR stations 102 and individually as a DGPSRstation 102) correct GPS errors to achieve precision GPS measurements atlocal rovers or vehicles 104 (referred to collectively asrovers/vehicles 104 and individually as a rover/vehicle 104). Accuracyof the GPS measurements may be on the order of several centimeters (cm).Each DGPSR station 102 can be positioned at a precisely known (e.g.,surveyed) location. The DGPSR station 102 measures its location asdetermined by GPS satellite signals from GPS satellites 106. Each DGPSRstation 102 computes the differential time error from its known locationrelative to its GPS-measured location. This differential error isconverted to equivalent differential range corrections, or pseudorangecorrections (PRC), for each of the satellites 106 in view. Thesedifferential corrections are sent via another radio channel 108 to thelocal vehicle/rover 104 to correct position calculations performed by aGPS receiver carried on the rover. The error drift rate of the PRC, orRange Rate Correction (RRC), may also be conveyed so the PRC can beappropriately corrected for a short time after the PRC is measured. Theterm “rover” or “rover platform” is construed broadly to include anyvehicle, platform, or device, whether fixed or mobile, whetherland-based, airborne, or seaborne. A rover or rover platform mayinclude, but is not limited to, a car, truck, tractor, train, airplane,ship, autonomous vehicle, fixed structure (e.g., tower or building),computer system, and the like.

The DGPSR station corrections require frequent updates (e.g., at leastonce a minute) due to changes in the differential error values overtime. Higher update rates (e.g., once per second) and other techniquesmay be needed to maintain greater RTK accuracy. The accuracy alsodegrades with distance from the DGPSR station location, maintaining highaccuracy over a limited range (e.g., up to 10 km). The local DGPSRstation 102 is useful for applications such as farming where autonomousfarm vehicles can operate over a limited area. Temporary DGPSR stationsare also useful for surveying.

Networked DGPS

FIG. 2 is an illustration of an environment for DGPS applicationscovering a wide area (e.g., North America). The environment of FIG. 2uses a network subscription service orchestrated at a DGPS serviceprovider 202. In FIG. 2, DGPSR stations 102 spanning the world arenetworked to feed differential corrections 203 into an Aggregated DGPS(ADGPS) database (DB) 204 that offers wide-area or even global DGPScorrection data. Since the DGPSR stations 102 are not sufficientlyclose-spaced over most of the potential coverage areas, interpolationand other methods are used to cover most regions. These evolvingservices might be offered to users (e.g., rovers 104) for a subscriptionfee for each device.

DGPS service provider 202 retrieves from ADGPS database 204 global DGPScorrection data, including clock and ephemeris errors, and differentialcorrections 203. DGPGS service provider 202 conveys, via 2-waycommunication links 206, differential corrections derived from theglobal DGPS correction data retrieved from the ADGPS database 204 to awireless modem and DGPS processor connected to a GPS receiver on eachrover 104. The rover 104 maintains 2-way contact with the DGPS serviceprovider 202 via the 2-way communication links to obtain frequent DGPSupdates (i.e., differential correction updates). For example, the rover104 first sends to the DGPS service provider 202 its (i.e., the rover's)approximate location as measured by the rover's GPS receiver. The DGPSservice provider 202 then computes the differential correctionsspecifically for the rover's location, and then sends the computeddifferential corrections to the rover. The differential correctionsinclude differential corrections for each of the satellites 106potentially in view at the rover's location. This may include up toabout 20 satellites for both L1 and L2 GPS frequencies updated at leastonce a minute.

The DGPS service provider 202 may compute the correction data (includingdifferential corrections) for the rover 104 by first establishing asingle virtual DGPSR station at the rover's reported location. Incontrast to an installed DGPSR station, a virtual DGPSR station is animaginary, unoccupied reference station having a virtual location, andfor which (virtual) differential correction data is created fromexisting differential correction data from installed DGPSR stationssurrounding the virtual location, as though the (virtual) differentialcorrection data had been observed by an installed DGPSR station at thevirtual location. The differential correction from the single virtualDGPSR station can be computed by interpolating the DGPS correction datafrom the closest (installed) DGPSR stations 102. The DGPS serviceprovider 202 may also use additional information about each satellite'sclock and ephemeris errors, as well as ionosphere and tropospherepropagation delays, to improve the correction information unique to therover's approximate location. Position accuracy of better than 10 cm hasbeen demonstrated using this technique.

A significant drawback of the networked DGPS services described inconnection with FIG. 2 is that the services require continuous, repeated2-way communications with the rover 104. Such 2-way communications maynot be available at all times and locations (coverage), and the user ofthe rover is likely to incur a subscription cost (either directly orindirectly). There is also a potential privacy issue since the serviceprovider knows the subscriber rover's location at all times.

Broadcasting a Differential Correction Map (DCM) to Rovers

Embodiments presented below, referred to as Differential Correction Map(DCM) embodiments, mitigate the drawbacks of the networked DGPSsubscription service as mentioned above. The DCM embodiments arepresented mainly in the context of GPS; however, the DCM embodimentsapply generally to, and work equally well in, any GNSS. The DCMembodiments include broadcasting a DCM to rovers by radio over a radiobroadcast coverage area (also referred to simply as a broadcast area),instead of using the limited, single-location virtual DGPS correctiondescribed above, which obviates the need for the constant 2-waycommunication between the rover and the DGPS service 202. The DCMrepresents a compact two-dimensional/surface description of the DGPScorrection data over/across the broadcast area. The surface may beplanar or the surface may be non-planar. An efficient broadcastmechanism includes broadcast radio, e.g., HD Radio. However, thebroadcast could also be over satellite, cellular network, or Internetcloud such as DTS Connected Radio. The broadcasting makes the DCMavailable to all rovers equipped with a radio receiver to receive thebroadcast that are located in the broadcast area. Consequently, the2-way service is replaced or augmented with the more efficient one-waybroadcast of the DCM covering the rover's location.

FIG. 3A is an illustration of an example environment that supportsbroadcasting of DGPS corrections to rovers via a DCM. The environment ofFIG. 3A includes the DGPSR stations 102, the rovers 104, and the ADGPSdatabase 204 mentioned above. Additionally, the environment of FIG. 3Aincludes a DCM DGPS (DDGPS) management entity or service 302 to generatea DCM and assist with disseminating the DCM, a radio station 304 totransmit a radio signal that conveys the DCM across a known broadcastarea 305 encompassing the rovers 104, and a communication network 306connected to the DDGPS management entity, the ADGPS database, and theradio station, and over which the aforementioned entities maycommunicate. DDGPS management entity 302 (also referred to more simplyas management entity 302) may include a service application hosted onone or more computer devices, such as servers, in a cloud-basedenvironment, to implement DCM methods presented herein. Communicationnetwork 306 may include wide area networks, such as the Internet, andlocal area networks.

An example high-level DCM broadcast method is now described briefly.Management entity 302 generates a constellation of virtual DGPSreference stations (not shown in FIG. 3A) having virtual locationsspaced-apart across the broadcast area of radio station 304. Moregenerally, the constellation may include/identify the virtual DGPSreference stations, or physical DGPS reference stations, or both. Theconstellation may include any number of such reference stations; in anexample, the constellation includes 4 or more reference stations.Management entity 302 stores information defining the constellation(i.e., stores the constellation) in a constellation database CDB hostedon a server, for example, along with other previously generatedconstellations.

Management entity 302 retrieves the above-mentioned global DGPScorrection data from the ADGPS database 204 over network 306. Managemententity 302 computes differential corrections for the virtual DGPSRstations of the constellation based on the information retrieved fromthe ADGPS database 204 and the virtual locations, and then generates theDCM based on the differential corrections for the constellation.Management entity 302 stores the DCM in a DCM database DDB hosted on theserver, for example, along with other previously computed DCMs.Management entity 302 sends the DCM to radio station 304 overcommunication network 306. In turn, radio station 304 transmits a radiosignal carrying the DCM across the broadcast area 305. While depicted asseparated, individual entities in FIG. 3A, it is understood thatmanagement entity 302 and radio station 304 may, in some embodiments, beimplemented as a single, combined entity at a given location.

Rover 104 includes a radio receiver configured to receive the radiosignal carrying the DCM, and to recover the DCM from the radio signal.The rover 104 includes a GPS receiver to receive satellite signals anddetermine an approximate location of the receiver (i.e., the rover)based on the satellite signals. The rover 104 evaluates the DCM at theapproximate location, to produce DCM-based differential corrections forthe satellite signals at the approximate location. The rover 104corrects the satellite signals using the DCM-based differentialcorrections, and determines a more accurate location for the rover basedon the (DCM) corrected satellite signals. For example, the roverincludes a DGPS processor to evaluate the received DCM at its (i.e., therover's) own GPS-measured location to determine its unique differentialcorrection data. This differential correction data can be converted tostandard RTCM format to conveniently replace the differential correctiondata delivered to typical RTK GPS receivers.

A single radio station is shown in FIG. 3A by way of example, only;however, it is understood that embodiments may include manygeographically separated radio stations connected to the managemententity 302 to support dissemination of geographically relevant DCMsacross geographically separated broadcast areas, as shown in FIG. 3B.FIG. 3B is an illustration of an example environment that supportsbroadcasting differential corrections to rovers across multiplebroadcast areas via multiple DCMs. The environment of FIG. 3B includesmultiple geographically separated radio stations 310(1)-310(N) totransmit respective radio signals in geographically separated respectivebroadcast areas 312(1)-312(N) encompassing respective rovers (notshown). Management entity 302 generates respective DCMs DCM(1)-DCM(N)for respective ones of broadcast areas 312(1)-312(N) based on respectiveconstellations of virtual DGPSR stations determined for/in each of thebroadcast areas. In an example, management entity 302 may selectexisting constellations from constellation database CDB, and generatethe DCMs based on corresponding ones of the selected DCMs. Also,management entity may select the existing DCMs from DCM database DDB.Management entity 302 sends DCM(1)-DCM(N) to respective ones of radiostations 310(1)-310(N), which in turn transmit the DCMs acrossrespective broadcast areas 312(1)-312(N).

FIG. 3C is an illustration of an example environment in which DCMs aremade accessible over communication network 306. In the environment ofFIG. 3C, rover 104 establishes a wireless network data connection tonetwork 306 based on any known or hereafter developed wireless standard,e.g., the Institute of Electrical and Electronics Engineers (IEEE)802.11 standards (i.e., WiFi) or the cellular standards, such as theLong Term Evolution (LTE) or other cellular standards. Then, managemententity 302 provides to rover 104 data packets including one or more DCMsover network 306 and the wireless connection between the rover and thecommunication network. Management entity 302 provides to rover 104 theDCM that covers the geographical area in which the rover is located.Management entity 302 and rover 104 may communicate with each using anyknown or hereafter developed communication protocol, including but notlimited to TCP/IP. In one example, management entity 302 may send theDCM to rover 104 in response to a request from the rover for the DCM. Inanother embodiment, management entity 302 may send the DCM to rover 104when the management entity detects, through signalling messagesoriginated in network 306, (i) that the rover has connected to thenetwork, and (ii) a general location of the rover indicated in thesignaling.

Generating the DCM and providing an efficient mechanism for representingthe DCM are described next.

DCM Generation

Management entity 302 can generate the DCM using differential correctiondata from multiple (i.e., a constellation of) virtual DGPSR stationsappropriately spaced over the DCM broadcast area, i.e., the broadcastarea of a given radio station. Since actual DGPSR stations sparselycover many areas, virtual DGPSR stations can be generatedconveniently/mathematically from the same DGPS database used bynetworked DGPS service providers, e.g., the ADGPS database 204. Thebroadcast area of a typical frequency modulation (FM) broadcast stationcan span up to about a 100 km radius. Variations in the differentialcorrections are assumed to be smooth over this area, so the differentialcorrections can be represented using a somewhat small number of(virtual) DGPSR stations spanning the DCM broadcast area. Presentmethods use a single virtual DGPSR station to provide a localcorrection, although it is not sufficient to provide an accuratedifferential correction over the broadcast area. Three virtual DGPSRstations define a differential correction plane over the broadcast areawhich can be useful if the change in differential varies approximatelylinearly with a change in distance, although generally more than threeDGPSRs are needed to accurately cover the typical broadcast area.

Several possibilities exist for defining the virtual DGPSR stations. Inthe embodiments presented herein, DGPSR stations can be located nearareas of greater use by rovers and roads with high rover density.Alternatively, virtual DGPSR stations can be arranged in one of a smallset of predetermined constellations. For example, the DCM could bedescribed using the differential correction data from 6 virtual DGPSRstations having virtual locations spaced over the broadcast area havinga radius of 100 km. Any geographic arrangement of virtual DGPSRs over abroadcast area is called a Virtual Reference Station Constellation(VRSC).

FIG. 4 is an illustration of an example VRSC 402 including 6 virtualDGPSR stations V0-V5 located at corresponding virtual locations, wherethe X and Y axes are in km. An HD radio station/transmitter is assumedto be near the virtual location of center virtual DGPSR station V0. Theother 5 virtual DGPSR stations V1-V5 are equally spaced from the centerby 100 km, which may be larger or smaller depending on a larger orsmaller broadcast area.

The data used to represent a first type of DCM referred to as a DCM Type0 include differential correction values (dT, a proxy for PRC) for eachfrequency and each satellite in view, along with the satelliteidentification for each of the DGPSR stations. Additionally, thelocations of the reference stations are conveyed to the rover with theDCM. However, if a common known VRSC is used, then the rover knows therelative locations within that known constellation, and only theabsolute position of the constellation (e.g., center virtual DGPSRstation location V0) is conveyed to the rover.

The differential corrections should be updated at least several timesper minute. For the example of FIG. 4, having 6 virtual DGPSRs in acommon known pattern and 20 satellites with frequencies L1 and L2, theminimum DCM size can be estimated at roughly 500 bytes, assuming 2 bytesfor each value. Assuming the DCM is broadcast a few times per minute,then the minimum throughput is on the order of a few hundred bits persecond.

Since the differential corrections can drift over time, the periodbetween DCM updates should be less than a minute. However, DCMinterruptions due to broadcast reception outages could result inexcessive errors in the differential corrections. This can be mitigatedby broadcasting the drift rate (dT, a proxy for RRC) of the differentialcorrections along with the DCM. This drift can be estimated by the firstderivative, or last difference, of the correction values. Each of thesevalues is small and requires less than one byte.

FIG. 9 is an illustration of an example DCM Type 0 (also referred to asa “Type 0 DCM”) that includes an identifier (ID) of the DCM (e.g., 0001,0002, and so on), a GNSS time (e.g., GPS time), and a location (e.g.,latitude/longitude) of a center of a reference station constellation,e.g., center location of reference station V0 of the VRSC describedabove. The DCM Type 0 may also convey locations of other referencesstations, e.g., of other virtual DGPSRs of the VRSC, either byidentifying a common known reference station pattern or the individualreference station locations, e.g., the virtual DGPSR station locationsfor reference stations V1-V5 for the VRSC pattern. In addition, the DCMType 0 also includes the differential correction values for eachreference station (e.g., each of reference stations V0-V5), eachfrequency (e.g., L1, L2), and each of the GPS satellites (e.g.,satellite IDs 0-19) in view, as well as some additional information. Inan example, each differential correction value may be an equivalent timedifferential dT represented in centimeters.

Proximity-Weighted Interpolation

As mentioned above, the rover (e.g., rover 104) uses DCM Type 0 tocompute the DGPS corrections (also referred to as DCM-based differentialcorrections) for its approximate location as given by the rover's GPSreceiver. In other words, the rover evaluates the DCM at the approximatelocation, to produce the DCM-based DGPS corrections. It is assumed thatthe DCM is comprised of the data necessary to describe the differentialerrors measured by a set of virtual DGPSR stations that span thebroadcast area, as described in FIG. 9. Instead of simply usingdifferential corrections determined by the single nearest virtual DGPSRstation location, the rover uses an interpolation method in which thedifferential correction at the approximate location of the rover isinfluenced by the respective differential corrections of all of thenearby virtual DGPSR stations.

To perform the interpolation, the rover computes a proximity-weightedcombination of the respective differential corrections (values) providedin DCM Type 0 influenced by the distance of the approximate locationfrom each of the virtual DGPSR stations. Mathematically, the rover'sinterpolated differential correction diff_(s) for satellite s can becalculated using the K differential correction values dT_(s,k) for eachof S satellites in view from the K virtual DGPSR stations identified inthe DCM. The rover's interpolated differential corrections for each ofthe S satellites can be computed as

${{{diff}_{s} = {\sum\limits_{k = 0}^{K - 1}{w_{k} \cdot {dT}_{s,k}}}};\mspace{25mu} {s = 0}},1,{{\ldots \mspace{14mu} S} - 1}$

Define the location of the k^(th) virtual DGPSR station having longitudeand latitude coordinates (x_(k), y_(k)). Next define the weight Wk usedfor the interpolation of the k^(th) virtual DGPSR station's differentialcorrection. This weight is a function of the reciprocal of its distancefrom the rover location (x,y). A useful expression for the weight Wk isrelated to the distance from the rover to the k^(th) virtual location(x_(k),y_(k)).

w_(k) = v_(k) ⋅ Norm${{{where}\mspace{14mu} v_{k}} = \frac{1}{\left\lbrack {\left( {x - x_{k}} \right)^{2} + \left( {y - y_{k}} \right)^{2}} \right\rbrack^{p} + ɛ}};{and}$${Norm} = \frac{1}{\sum\limits_{k}^{K - 1}v_{k}}$

Each weight wk is a function of the distance squared, to the power p,plus E. The power p is used to control the sensitivity of theinterpolation to the nearest virtual DGPSR station; the recommendedvalues for p are 1 or 2. The addition of e prevents overflow due todivision by zero. The intermediate weights v_(k) are normalized by Normsuch that the weights wk sum to one, preventing a bias factor for dT.

FIGS. 5A and 5B are example isometric DCM plots diff_1 and diff_2 for 6interpolated DGPSR differential corrections for a given satellite takenfrom the example of FIG. 4, for p=1 and p=2, respectively. Arbitrarydifferential corrections from the virtual DGPSR stations are chosen forthis example. The center (V0) DGPSR differential value is set todT_(s,0)=100, and the other five locations to dT_(s,k)=90, representingboth plots. The units of time for values of dT and diff in this examplecould be equivalent to centimeters of propagation distance (range) atthe speed of light.

Management entity 302 can improve overall DCM accuracy by establishing agreater number of virtual DGPSR stations in the constellation. Noticethat the plots of FIGS. 5A and 5B show some sensitivity or distortionbetween the locations of the virtual DGPSR stations. It is not reallyknown what the differential correction values should be without morevirtual locations or additional information. However, additional virtuallocations increase the amount of DCM data and required transmissionthroughput. One way to maintain the DCM size while increasing accuracyis to alternate VRSCs and corresponding DCMs on successive transmissionsof the DCMs. That is, management entity 302 generates multiplealternating constellations covering the broadcast area (andcorresponding DCMs), as shown in FIG. 6A. FIG. 6A shows the virtualDGPSR stations of a first constellation (e.g., the constellation 402 ofFIG. 4) as separated circles, and the virtual DGPSR stations of asecond/alternate constellation as squares. Notice that the alternatevirtual locations (i.e., squares) are spaced between the originalvirtual locations (i.e., circles) around the outer perimeter, and thecenter virtual location is repeated. This pair of constellations hascomplementary symmetric coverage. The alternating constellations resultin alternating DCMs, one per constellation, for transmission.

With respect to alternating DCM constellations, upon startup, the roveruses the first DCM constellation it receives. The rover can averageinterpolation results from subsequent alternating DCMs for even betterresolution accuracy. More generally, this could apply to severalalternate constellations, and lossy averaging with an Infinite ImpulseResponse (IIR) filter could be used. Management entity 302 can identifyeach constellation of the set of alternate constellations with aconstellation/DCM identifier of one or more bits in a header/controlfield of the DCM, assuming a standard constellation. FIG. 6B is anexample isometric DCM plot diff_avg resulting from averaging theconsecutive complementary constellations shown in FIG. 6A at the rover.Notice the apparent reduced sensitivity or distortion between thereference locations.

Polynomial Interpolation

The Proximity-weighted interpolation described above in connection withconstellation 402 of FIG. 4 includes the known locations of the virtualDGPSR stations V0-V5 and their differential errors. The roverinterpolates between the differential errors of DCM Type 0.

Alternatively, the previously described DCM Type 0 of FIG. 9 used inproximity-weighted interpolation can be replaced with a second type ofDCM, referred to as a DCM Type 1 (also referred to as a “Type 1 DCM”),which is a polynomial-fit DCM, described below in connection with FIG.10. DCM Type 1 could use the same 6 virtual locations discussed above,or different ones. DCM Type 1 for satellite s of the previous exampleusing 6 virtual DGPSR stations, V0 through V5, can be described using 6scalar polynomial coefficients A, B, C, D, E, F for each frequency ofeach satellite in view.

diff_(s) =A _(s) +B _(s) ·x+C _(s) ·y+D _(s) ·x ² +E _(s) ·y ² +F _(s)x·y

The polynomial coefficients can be derived using the virtual DGPSRlocations and differential corrections dT. The 6 (x_(k),y_(k))coordinates form a linear system of 6 two-dimensional quadraticequations to solve for the polynomial coefficients for each satellite inview.

A+B·x ₀ +C·y ₀ +D·x ₀ ² +E·y ₀ ² +F·x ₀ ·y ₀ =dT ₀

A+B·x ₁ +C·y ₁ +D·x ₁ ² +E·y ₁ ² +F·x ₁ ·y ₁ =dT ₁

A+B·x ₂ +C·y ₂ +D·x ₂ ² +E·y ₂ ² +F·x ₂ ·y ₂ =dT ₂

A+B·x ₃ +C·y ₃ +D·x ₃ ² +E·y ₃ ² +F·x ₃ ·y ₃ =dT ₃

A+B·x ₄ +C·y ₄ +D·x ₄ ² +E·y ₄ ² +F·x ₄ ·y ₄ =dT ₄

A+B·x ₅ +C·y ₅ +D·x ₅ ² +E·y ₅ ² +F·x ₅ ·y ₅ =dT ₅

which can be more-compactly represented in matrix form as

${Q \cdot \begin{pmatrix}A \\B \\C \\D \\E \\F\end{pmatrix}} = {{{dT}\mspace{14mu} {Where}{\mspace{14mu} \;}Q} = {{\begin{pmatrix}1 & x_{0} & y_{0} & x_{0}^{2} & y_{0}^{2} & {x_{0} \cdot y_{0}} \\1 & x_{1} & y_{1} & x_{1}^{2} & y_{1}^{2} & {x_{1} \cdot y_{1}} \\1 & x_{2} & y_{2} & x_{2}^{2} & y_{2}^{2} & {x_{2} \cdot y_{2}} \\1 & x_{3} & y_{3} & x_{3}^{2} & y_{3}^{2} & {x_{3} \cdot y_{3}} \\1 & x_{4} & y_{4} & x_{4}^{2} & y_{4}^{2} & {x_{4} \cdot y_{4}} \\1 & x_{5} & y_{5} & x_{5}^{2} & y_{5}^{2} & {x_{5} \cdot y_{5}}\end{pmatrix}{\mspace{11mu} \mspace{14mu}}{and}\mspace{14mu} {dT}} = {{\begin{pmatrix}{dT}_{0} \\{dT}_{1} \\{dT}_{2} \\{dT}_{3} \\{dT}_{4} \\{dT}_{5}\end{pmatrix}\mspace{14mu} {Then}\mspace{14mu} \begin{pmatrix}A \\B \\C \\D \\E \\F\end{pmatrix}} = {Q^{- 1} \cdot {dT}}}}}$

Q is a 6 by 6 matrix of elements derived from the VRSC coordinates, V0through V5, in the DCM polynomial.

FIG. 7A presents an isometric DCM plot diff_poly for the constellation402 of FIG. 4 using a two-dimensional quadratic polynomial fit (alsoreferred to simply as a “polynomial”) representation of DCM Type 1. Theradio receiver (and associated DCM processing) at the rover needs noknowledge of the virtual DGPSR locations using polynomialrepresentation, other than the absolute location of the VRSC (e.g., V0).So any convenient VRSC can be used in creating the DCM polynomial.Although DCM Type 1 passes through all the reference points, the distantportions of the DCM tend to increase or decrease “quadratically” withthe square of the distance from center. Within the broadcast area, thismay be mitigated by limiting the computed radial distance from thecenter location of the DCM. This modification is shown in the isometricDCM plot diff_poly_mod of FIG. 7B, which shows that the isometric valuesare lower bounded to about 90 cm (the lowest of the DGPSR values).

The polynomial representation of DCM Type 1 can be converted to anequivalent DCM Type 0 representation of the VRSC shown in FIG. 9. Thisis accomplished simply by evaluating the polynomial at the virtual DGPSRlocations. Conversely, it was already shown that a polynomialrepresentation can be derived from a VRSC. However, at least K virtualDGPSR stations are preferably used to derive a polynomial with Kcoefficients. Furthermore, the useful choices of K for 2-dimensionalpolynomials are 3 (linear), 6 (quadratic), 10 (cubic), 15 (quartic), 21(quintic), etc.

Both the proximity-weighted and the polynomial interpolation methodsproduce different interpretations of the DCM. Although the 2interpolation methods produce DCMs having the same correction values(diff) at the virtual DGPSR locations, the resulting DCMs differ atother locations. Overall DCM accuracy can be improved by adding morevirtual locations or higher-order polynomials. The accuracy of DCM Type0 interpolation using the proximity-weighted method can be improved byadjusting the power p in the interpolation expression. Additionalreference locations are then used to determine which power p offers abetter fit. This can be attractive since the only additional informationconveyed by the DCM is p.

The accuracy of the polynomial representation of DCM Type 1 can beimproved without increasing its order above K by using least mean square(LMS) methods to fit more than K reference locations, although thereference locations will no longer be exact solutions to the polynomial.Furthermore, overall error can be reduced by averaging polynomials ofalternate VRSCs with interlaced location spacing at DCM creation. Noadditional information needs to be conveyed with the DCM in this case.

The DCM accuracy can also be improved by using a weighted average ofboth the proximity-weighted and polynomial interpolation methods. Anadditional weighting parameter a that reduces the overall DCM error canbe conveyed with the DCM to enable this averaging at the receiver. Aweighted average example (α=0.5) for the single 6-location constellationexample is shown in FIG. 8. In other words, FIG. 8 is an illustration ofan example isometric DCM plot for an average fit of the modifiedpolynomial and proximity-weighted interpolation methods.

diff_(s)=α·diff_prox_(s)+(1−α)−diff_poly_(s)

Encryption of DCM

Broadcast of the DCM makes the service available to all rovers withinthe broadcast area. The DCM can be encrypted if conditional access isdesired; for example if a subscription service is offered. Encryptioncan be implemented by assigning each receiver a unique code (e.g.,serial number). A monthly encryption key could be used to updatesubscribers through internet access.

Flowcharts

FIGS. 11, 12, and 13A-13D present flowcharts of various methods directedto DCMs drawn from operations described above.

FIG. 11 is a flowchart of an example method 1100 of generating anddisseminating a DCM in any of the environments of FIGS. 3A-3C.

At 1102, management entity 302 determines/establishes a referencestation constellation of GNSS reference stations for a geographicalarea. To do this, management entity 302 may initially define/generatethe reference station constellation, or the management entity may selectthe reference station constellation from among a set of predeterminedreference station constellations. If the determined reference stationconstellation is not presently stored in memory, management entity 302stores the reference station constellation to memory, e.g., toconstellation database CDB. The reference stations of the referencestation constellation may include only physical reference stations, onlyvirtual reference stations, or a mix of both. In one arrangement, thereference stations have locations that are non-uniformly spaced acrossthe geographical area. In another arrangement, the references stationshave locations that are approximately uniformly spaced across thegeographical area. In the latter arrangement, the reference stations mayall be virtual reference stations. In an embodiment, the referencestation constellation includes at least four reference stations.

At 1104, management entity 302 access/retrieves from ADGPS database 204aggregated differential correction information corresponding to thegeographical area.

At 1106, management entity 302 generates a DCM representing a surface ofdifferential corrections (e.g., dTs) for the GNSS, wherein the surfaceof the differential corrections varies across the geographical arearepresented by the DCM. In one embodiment, the surface of thedifferential corrections may be planar. In another embodiment, thesurface of the differential corrections may be non-planar in that afirst derivative across any line (e.g., latitude or longitude) over thesurface may not be constant, e.g., is not constant. The differentialcorrections are based on the reference station constellation ofreference stations (which may include the physical or the virtualreference stations) having respective locations spanning thegeographical area, and the information accessed from ADGPS database 204.The surface of the DCM is configured to be evaluated at any location inthe geographical area, to produce DCM-based differential corrections forthat location.

Management entity 302 may generate the DCM as a DCM Type 0, which is arepresentation of the reference station constellation. The DCM Type 0includes a DCM identifier (ID), GNSS time, the locations of one or moreof the reference stations, and differential corrections for eachreference station. More specifically, the DCM includes information asshown in FIG. 9, such as the GNSS time, locations of the referencestations, an identifier of each GNSS satellite in view of thegeographical area represented by the DCM, and the differentialcorrection for each reference station for each frequency of each GNSSsatellite in view. The DCM may further include a drift rate of thedifferential correction (e.g., dT) for each reference station for eachfrequency of each GNSS satellite in view. This drift rate can beestimated by the first derivative (dT), or last difference, of thedifferential correction values.

Management entity 302 may generate the DCM as a DCM Type 1, in whichcase the management entity generates the DCM from a two-dimensionalpolynomial surface fit to the differential correction associated witheach frequency and each GNSS satellite in view of the geographical arearepresented by the DCM. The DCM Type 1 is an indirect representation ofthe reference station constellation compared to the DCM Type 0. The DCMType 1 includes a DCM identifier (ID), GNSS time, the locations of oneor more of the reference stations, and polynomial coefficients of thepolynomial surface fit for the differential corrections for thereference stations. More specifically, the DCM Type 1 includesinformation as shown in FIG. 10, which shows an example DCM Type 1 basedon the polynomial fit corresponding to the constellation 402 of FIG. 4.The DCM Type 1 includes the GNSS time, an origin of the two-dimensionalpolynomial surface fit (e.g., the location of reference station V0), anidentifier of each of the GNSS satellites (e.g., satellite IDs SAT ID0-19) in view of the geographical area represented by the DCM,polynomial coefficients for the reference stations associated with thetwo-dimensional polynomial surface fit (e.g., polynomial coefficientsA-F for reference stations V0-V6, respectively). The DCM may alsoinclude a drift rate of each polynomial coefficient for each frequencyof each GNSS satellite in view. This drift rate can be estimated by thefirst derivative, or last difference, of the polynomial coefficients

Management entity optionally encrypts the DCM using any known orhereafter developed encryption technique.

At 1108, management entity 302 makes the DCM (in encrypted orunencrypted form) accessible to interested entities. In one embodiment,management entity 302 makes the DCM accessible over communicationnetwork 306, in which case a rover 104 may connect to the network andaccess the DCM over the network. In another embodiment, managemententity 302 provides the DCM to a radio transmitter (e.g., radiotransmitter 304) that broadcasts the DCM across the geographical areavia a radio signal, i.e., the radio signal is modulated to convey theDCM. This enables rovers 104 in the geographical area to receive the DCMvia the radio signal. Example broadcast radio outlets for transmittingthe DCM include, but are not limited to, FM subsidiary communicationsauthorization (SCA) digital subcarrier (e.g., RBDS), digital radiobroadcast (e.g., HD Radio, digital audio broadcast (DAB), Digital RadioMondiale (DRM), and the like), digital satellite broadcast (e.g., SiriusXM), and digital TV broadcast.

The information in ADGPS database 204 and/or the constellation used in agiven coverage are dynamic and thus may change periodically.Accordingly, method 1100 repeats over time to periodically update theDCM so that the DCM accurately reflects such changes.

When rovers 104 located in the geographical area receive the DCM, therovers may evaluate the received DCM at their respective locations toderive accurate locations for the rovers, as described herein. Forexample, a given rover may use a GPS receiver to derive an approximatelocation of the rover based on GPS satellite signals, evaluate the(received) DCM at the approximate location to produce DCM-baseddifferential corrections at the approximate location, and then correctthe satellite signals using the DCM-based differential corrections, toproduce a more accurate location of the rover.

Variations of method 1100 are possible. In one variation, at operation1102, management entity 302 selects a VRSC from a predetermined set ofVRSCs, and generates the DCM from the selected predetermined VRSC.

In another variation, management entity 302 generates multiple DCMsbased on multiple reference station constellations (i.e., one DCM percorresponding VRSC) for approximately the same geographical area, andalternately transmits the multiple DCMs. For example, management entitytransmits a first DCM, a second DCM, a third DCM, and then the firstDCM, the second DCM, the third DCM, and so on.

In yet another variation, management entity 302 generates multiple DCMsincluding a DCM Type 0 and a DCM Type 1 for the same geographical area,and then transmits both DCMs, e.g., alternately transmits the DCMs.

Generate DCM operation 1106 is described above in connection with otheroperations of method 1100 by way of example, only. It is understood thatgenerate DCM operation 1106, taken alone, represents a stand-aloneoperation or method. In other words, operation 1106 alone represents itsown single-operation method for which operations 1102, 1104, and 1108are unnecessary or not essential.

FIG. 12 is a flowchart of an example method 1200 of generating andbroadcasting multiple DCMs across multiple broadcast areas performed inthe network environment of FIG. 3B.

At 1202, management entity 302 generates multiple reference stationconstellations, each corresponding to a respective one of broadcastareas 312(1)-312(N). Each constellation includes reference stationslocated in the broadcast area to which the constellation corresponds.

At 1204, management entity 302 generates multiple distinct DCMsDCM(1)-DCM(N), each uniquely associated with a respective one ofmultiple broadcast areas 312(1)-312(N), and each based on a respectiveone of the multiple reference station constellations. Management entity302 provides the DCMs DCM(1)-DCM(N) to respective ones of radio stations310(1)-310(N), which transmit the DCMs into their respective broadcastareas.

FIG. 13A is a flowchart of an example method 1300 of processing the DCM,performed by a processor (e.g., processor 1506) of a rover platform(e.g., rover 104). The rover platform may also include GPS receiver1502, radio receiver 1504, and/or wireless network interface unit (NIU)1505, as well as memory, coupled to the processor, as described below inconnection with FIG. 15.

At 1302, the processor receives an approximate location of the roverplatform based on satellite signals for a GNSS. For example, GPSreceiver 1502 may receive the satellite signals, determine theapproximate location of the rover platform based on the satellitesignals (e.g., time delays associated with the satellite signals), andprovide the approximate location to the processor.

At 1304, the processor receives a DCM. For example, radio receiver 1504receives a broadcast radio signal that conveys the DCM, and provides theDCM to the processor. Alternatively, NIU 1505 receives data packets thatconvey the DCM over a wireless network connection, and sends the DCM tothe processor. The DCM represents a non-planar surface of differentialcorrections for the GNSS that varies across a geographical arearepresented by the DCM. The differential corrections are based on areference station constellation of GNSS reference stations havingrespective locations spanning the geographical area. The GNSS referencestations may include physical and/or virtual reference stations.

A first type of DCM includes a DCM identifier (ID), GNSS time, thelocations of one or more of the reference stations, and differentialcorrections for each reference station.

A second type of DCM includes a DCM identifier (ID), GNSS time, thelocations of one or more of the reference stations, and polynomialcoefficients of a polynomial fit for the differential corrections forthe reference stations.

At 1306, the processor derives DCM-based differential corrections forthe satellite signals at the approximate location based on the DCM.

For the first type of DCM, the processor interpolates the DCM at theapproximate location to produce the DCM-based differential correctionsas interpolated differential corrections. The interpolating includesproximity-weighting the differential corrections of the referencestations as a function of distances between the approximate location andthe locations of the reference stations, and combining theproximity-weighted differential corrections to produce the interpolateddifferential corrections.

For the second type of DCM, the processor evaluates the polynomial-fitat the approximate location using the polynomial coefficients.

At 1308, the processor corrects the satellite signals using theDCM-based differential corrections, to produce corrected satellitesignals.

At 1310, the processor determines a location of the rover platform usingthe corrected satellite signals. The processor may use the location inany number of navigation applications hosted on the rover, includingdisplay of the location.

FIG. 13B is a flowchart of an example method 1320 of processing multipleDCMs that expands on or is alternative to operations 1304 and 1306described above, performed by the processor of the rover platform.

At 1322, the processor receives multiple DCMs including a first DCM ofthe first type mentioned above in method 1300 and a second DCM of thesecond type mentioned above.

At 1324, the processor:

-   -   a. Derives first differential corrections based on the first DCM        using the interpolating described above;    -   b. Derives second differential corrections based on the second        DCM using the evaluating the polynomial described above; and    -   c. Computes a weighted average of the first differential        corrections and the second differential corrections, to produce        the DCM-based differential corrections as weighted-averaged        differential corrections.

FIG. 13C is a flowchart of an example method 1330 of processing multipleDCMs that expands on or is alternative to operations 1304 and 1306described above, performed by the processor of the rover platform.

At 1332, the processor receives multiple DCMs, which may be of the firsttype or the second type.

At 1334, the processor derives differential corrections based on eachDCM, and computes a weighted average of the differential correctionsfrom the multiple DCMs, to produce the DCM-based differentialcorrections as weighted-averaged differential corrections.

FIG. 13D is a flowchart of an example method 1340 of processing a DCMthat expands on or is alternative to method 1300.

At 1342, the processor receives a DCM of the first type or the secondtype, and which includes the location of only one of the referencestations of the constellation.

At 1344, the processor derives locations of other reference stations ofthe constellation from the location of the only one reference stationand known relative locations of other ones of the reference stationsrelative to the location of the only one reference station that arestored in the memory of the rover as a priori information.

Once all of the locations of the reference stations are known from 1342and 1344, the processor derives the DCM-based differential correctionsas described above.

Management Entity, Rover, and Transmit Station Block Diagrams

With reference to FIG. 14, there is shown an example block diagram forDDGPS management entity 302. In the example, management entity 302includes a computer system, such as a server, having one or moreprocessors 1410, a network interface unit (NIU) 1412, and a memory 1414.Memory 1414 stores control software 1416 (referred as “control logic”),that when executed by the processor(s) 1410, causes the computer systemto perform the various operations described herein for management entity302.

The processor(s) 1410 may be a microprocessor or microcontroller (ormultiple instances of such components). The NIU 1412 enables managemententity 302 to communicate over wired connections or wirelessly with anetwork, such as network 306. NIU 1412 may include, for example, anEthernet card or other interface device having a connection port thatenables management entity 302 to communicate over the network via theconnection port. In a wireless embodiment, NIU 1412 includes a wirelesstransceiver and an antenna to transmit and receive wirelesscommunication signals to and from the network.

The memory 1414 may include read only memory (ROM), random access memory(RAM), magnetic disk storage media devices, optical storage mediadevices, flash memory devices, electrical, optical, or other physicallytangible (i.e., non-transitory) memory storage devices. Thus, ingeneral, the memory 1414 may comprise one or more tangible(non-transitory) computer readable storage media (e.g., memorydevice(s)) encoded with software or firmware that comprises computerexecutable instructions. For example, control software 1416 includeslogic to implement operations performed by the management entity 302.Thus, control software 1416 implements the various methods/operationsdescribed above, including methods 1100 and 1200. Memory 1414 alsostores data 1418 generated and used by control software 1416, includingDCMs (e.g., DCM database DDB), information for constellations of virtualand/or physical DGPSR stations (e.g., constellation database CDB), andaggregated correction information, and so on.

A user, such as a network administrator, may interact with managemententity 302 through a user device 1420 (also referred to as a “networkadministration device”) that connects by way of a network withmanagement entity 302. The user device 1420 may be a personal computer(laptop, desktop), tablet computer, SmartPhone, and the like, with userinput and output devices, such as a display, keyboard, mouse, and so on.Alternatively, the functionality and a display associated with userdevice 1420 may be provided local to or integrated with managemententity 302.

FIG. 15 is an example block diagram of an apparatus carried by/affixedto the rover 104 for performing embodiments presented herein. Rover 104carries/includes a GPS receiver 1502 to receive and process GPSsatellite signals, a radio receiver 1504 to receive and process a radiosignal transmitted by radio station(s) 304 and/or 310(1)-310(N), forexample, a wireless network interface unit 1505 to establish 2-waywireless connections to a communication network, a processor 1506, and amemory 1508 all coupled to, and able to communicate with, one another.Wireless network interface unit 1505 may include a Wi-Fi interfaceand/or a cellular interface for transmitting and receiving wirelesssignals. In some embodiments, the network interface unit may alsoinclude a wired interface, such as an Ethernet card including anEthernet port. Radio receiver 1504 and wireless network interface unit1505 are examples of wireless receivers to receive signals wirelessly.

GPS receiver 1502 determines an approximate location of the GPS receiverbased on time delays of the satellite signals measured at the GPSreceiver. GPS receiver 1502 provides the approximate location and thetime delays to processor 1506. In an embodiment, GPS receiver 1502 andvarious functions performed by processor(s) 1506 may be combined in aDGPS receiver of the rover. Radio receiver 1504 or network interfaceunit 1505 recovers the DCM from the radio signal or from thecommunication network and provides the DCM to processor 1506. Radioreceiver 1504 and network interface unit 1505 may each include abaseband processor (BBP) to assist with processing, e.g., demodulatingand decoding, radio signals. Alternatively, the baseband processor maybe included with processor(s) 1506, or be separate from theprocessor(s), radio receiver, or network interface unit.

Memory 1508 stores control software 1510 (referred as “control logic”),that when executed by the processor 1506, causes the rover to performthe various operations described herein. Memory 1508 may include readonly memory (ROM), random access memory (RAM), magnetic disk storagemedia devices, optical storage media devices, flash memory devices,electrical, optical, or other physically tangible (i.e., non-transitory)memory storage devices. Thus, in general, the memory 1508 may compriseone or more tangible (non-transitory) computer readable storage media(e.g., memory device(s)) encoded with software or firmware thatcomprises computer executable instructions. For example, controlsoftware 1510 includes logic to implement operations performed by therover. Thus, control software 1510 implements the variousmethods/operations described above, including the methods connected withFIGS. 13A-13D. Memory 1508 also stores data 1512 generated and used bycontrol software 1510, including DCMs, DCM identifiers, approximatelocations, time delays, corrected time delays, DCM-based differentialcorrections for the time delays, and so on.

FIG. 16 is a block diagram of an example radio station 1600representative of any of radio stations 304 and 310(1)-310(N). In theexample, radio station 1600 includes a processor 1610, an NIU 1612, amemory 1614, and a radio transmitter 1616 all coupled to, and thus ableto communicate with, one another. Memory 1614 stores control software1618 (referred as “control logic”), that when executed by the processor1610, causes the processor to control the radio station 1600 to performthe various operations described herein. The processor 1610 may be amicroprocessor or microcontroller (or multiple instances of suchcomponents). The NIU 1612 enables radio station 1600 to communicate overwired connections or wirelessly with a network, such as network 306. NIU1612 may include, for example, an Ethernet card or other interfacedevice having a connection port that enables the radio station 1600 tocommunicate over the network via the connection port. In a wirelessembodiment, NIU 1612 includes a wireless transceiver and an antenna totransmit and receive wireless communication signals to and from thenetwork.

The memory 1614 may include read only memory (ROM), random access memory(RAM), magnetic disk storage media devices, optical storage mediadevices, flash memory devices, electrical, optical, or other physicallytangible (i.e., non-transitory) memory storage devices. Thus, ingeneral, the memory 1614 may comprise one or more tangible(non-transitory) computer readable storage media (e.g., memorydevice(s)) encoded with software or firmware that comprises computerexecutable instructions. For example, control software 1618 includeslogic to implement operations performed by the radio station 1600, suchas receiving messages conveying a DCM from the NIU 1612, and sending theDCM to radio transmitter 1616 for transmission. Memory 1614 also storesdata 1620 generated and used by control software 1618, including DCMs,radio station content for transmission, and so on.

Radio transmitter 1616 receives data/information (e.g., a DCM andprogram content for a radio station) from processor 1610, modulates acarrier signal with the data/information, and transmits/broadcasts themodulated carrier signal as a radio signal or radio frequency across abroadcast area. Radio station 1616 may include a baseband processor (notshown) to perform baseband modulation, or the baseband processor may beintegrated with processor 1610.

In some embodiments, a DGPS management entity may incorporate/combinethe full functionality of both radio station 1600 and management entity302, such that the DGPS management entity also includes a radiotransmitter (e.g., radio transmitter 1616) for transmission of DCMs.

In summary, in one embodiment, a method is provided comprising:generating for a Global Navigation Satellite System (GNSS) adifferential correction map (DCM) representing a non-planar surface ofdifferential corrections that varies across a geographical arearepresented by the DCM, the differential corrections being based on areference station constellation of GNSS reference stations havingrespective locations spanning the geographical area, the GNSS referencestations including physical or virtual reference stations.

In another embodiment, a method is provided comprising: determining avirtual reference station constellation (VRSC) of virtual referencestations having respective locations spanning a geographical area;generating for a Global Navigation Satellite System (GNSS) adifferential correction map (DCM) based on the VRSC, the DCMrepresenting a non-planar surface of differential corrections thatvaries across a geographical area represented by the DCM; andtransmitting a radio signal that conveys the DCM across the geographicalarea.

In yet another embodiment, an apparatus is provided comprising: memory;and a processor coupled to the memory and configured to: generate for aGlobal Navigation Satellite System (GNSS) a differential correction map(DCM) representing a non-planar surface of differential corrections thatvaries across a geographical area represented by the DCM, thedifferential corrections being based on a reference stationconstellation of GNSS reference stations having respective locationsspanning the geographical area, the GNSS reference stations includingphysical or virtual reference stations.

In a further embodiment, a non-transitory computer readable medium isprovided. The computer readable medium is encoded with instructionsthat, when executed by a processor, cause the processor to generate fora Global Navigation Satellite System (GNSS) a differential correctionmap (DCM) representing a non-planar surface of differential correctionsthat varies across a geographical area represented by the DCM, thedifferential corrections being based on a reference stationconstellation of GNSS reference stations having respective locationsspanning the geographical area, the GNSS reference stations includingphysical or virtual reference stations

In another embodiment, a method is provided comprising: receiving anapproximate location of a rover platform based on satellite signals fora Global Navigation Satellite System (GNSS); receiving for the GNSS adifferential correction map (DCM) representing a non-planar surface ofdifferential corrections that varies across a geographical arearepresented by the DCM, the differential corrections being based on areference station constellation of GNSS reference stations havingrespective locations spanning the geographical area, the GNSS referencestations including physical or virtual reference stations; derivingDCM-based differential corrections for the satellite signals at theapproximate location based on the DCM; correcting the satellite signalsusing the DCM-based differential corrections; and determining a locationof the rover platform using the corrected satellite signals.

In yet another embodiment, an apparatus is provided comprising: memory;and a processor coupled to the memory and configured to: receive anapproximate location of a rover platform that is based on satellitesignals for a Global Navigation Satellite System (GNSS); receive for theGNSS a differential correction map (DCM) representing a non-planarsurface of differential corrections that varies across a geographicalarea represented by the DCM, the differential corrections being based ona reference station constellation of GNSS reference stations havingrespective locations spanning the geographical area, the GNSS referencestations including physical or virtual reference stations; deriveDCM-based differential corrections for the satellite signals at theapproximate location based on the DCM; correct the satellite signalsusing the DCM-based differential corrections; and determine a locationof the rover platform using the corrected satellite signals.

In an even further embodiment, a non-transitory computer readable mediumencoded with instructions is provided. The instructions, when executedby a processor, cause the processor to: receive an approximate locationof a rover platform that is based on satellite signals for a GlobalNavigation Satellite System (GNSS); receive for the GNSS a differentialcorrection map (DCM) representing a non-planar surface of differentialcorrections that varies across a geographical area represented by theDCM, the differential corrections being based on a reference stationconstellation of GNSS reference stations having respective locationsspanning the geographical area, the GNSS reference stations includingphysical or virtual reference stations, the DCM including a DCMidentifier (ID), GNSS time, the locations of one or more of thereference stations, and differential corrections for each referencestation; interpolate the DCM at the approximate location to produceinterpolated differential corrections for the satellite signals; correctthe satellite signals using the interpolated differential corrections;and determine a location of the rover platform using the correctedsatellite signals.

Although the techniques are illustrated and described herein as embodiedin one or more specific examples, it is nevertheless not intended to belimited to the details shown, since various modifications and structuralchanges may be made within the scope and range of equivalents of theclaims.

Each claim presented below represents a separate embodiment, andembodiments that combine different claims and/or different embodimentsare within the scope of the disclosure and will be apparent to those ofordinary skill in the art after reviewing this disclosure.

What is claimed is:
 1. A method comprising: generating for a GlobalNavigation Satellite System (GNSS) a differential correction map (DCM)representing a non-planar surface of differential corrections thatvaries across a geographical area represented by the DCM, thedifferential corrections being based on a reference stationconstellation of GNSS reference stations having respective locationsspanning the geographical area, the GNSS reference stations includingphysical or virtual reference stations.
 2. The method of claim 1,wherein the reference station constellation includes a virtual referencestation constellation (VRSC) including virtual reference stationsapproximately uniformly spaced across the geographical area representedby the DCM.
 3. The method of claim 2, further comprising: selecting theVRSC from a set of predetermined VRSCs, wherein the generating includesgenerating the DCM for the selected predetermined VRSC.
 4. The method ofclaim 3, wherein the DCM includes: GNSS time; and a location of only oneof the reference stations; an identifier of the VRSC; an identifier ofeach GNSS satellite in view of the geographical area represented by theDCM; and the differential correction for each reference station for eachfrequency of each GNSS satellite in view.
 5. The method of claim 1,wherein the reference stations are non-uniformly spaced across thegeographical area represented by the DCM.
 6. The method of claim 1,wherein the DCM includes an identifier of the DCM, GNSS time, thelocations of one or more of the reference stations, and differentialcorrections for each reference station.
 7. The method of claim 1,wherein the DCM includes: GNSS time; locations of the referencestations; an identifier of each GNSS satellite in view of thegeographical area represented by the DCM; and the differentialcorrection for each reference station for each frequency of each GNSSsatellite in view.
 8. The method of claim 7, wherein the DCM furtherincludes a drift rate of the differential correction for each referencestation for each frequency of each GNSS satellite in view.
 9. The methodof claim 1, wherein the DCM includes an identifier of the DCM, GNSStime, the locations of one or more of the reference stations, andpolynomial coefficients of a polynomial fit for the differentialcorrections for the reference stations.
 10. The method of claim 1,wherein the generating includes generating the DCM from atwo-dimensional polynomial surface fit to the differential correctionassociated with each frequency and each GNSS satellite in view of thegeographical area represented by the DCM.
 11. The method of claim 10,wherein the DCM includes: GNSS time; an origin of the two-dimensionalpolynomial surface; an identifier of each of the GNSS satellites in viewof the geographical area represented by the DCM; and polynomialcoefficients associated with the two-dimensional polynomial surface. 12.The method of claim 11, wherein the DCM further includes a drift rate ofeach polynomial coefficient for each frequency of each GNSS satellite inview.
 13. The method of claim 1, further comprising transmitting the DCMover a broadcast coverage area via a radio signal.
 14. The method ofclaim 1, further comprising storing the DCM in a computer databaseaccessible over a communication network.
 15. The method of claim 1,further comprising: prior to the generating, accessing from a databaseaggregated differential correction information, wherein the generatingincludes generating the DCM from the aggregated differential correctioninformation.
 16. The method of claim 1, further comprising: generatingmultiple DCMs based on multiple reference station constellations, andalternately transmitting the multiple DCMs.
 17. The method of claim 1,further comprising: encrypting the DCM; and transmitting the encryptedDCM via broadcast radio signal or over a communication network.
 18. Amethod comprising: determining a virtual reference station constellation(VRSC) of virtual reference stations having respective locationsspanning a geographical area; generating for a Global NavigationSatellite System (GNSS) a differential correction map (DCM) based on theVRSC, the DCM representing a non-planar surface of differentialcorrections that varies across a geographical area represented by theDCM; and transmitting a radio signal that conveys the DCM across thegeographical area.
 19. The method of claim 18, wherein the VRSC includesfour or more of the virtual reference stations.
 20. The method of claim19, wherein the virtual reference stations are approximately uniformlyspaced across the geographical area represented by the DCM.
 21. Themethod of claim 19, wherein the virtual reference stations arenon-uniformly spaced across the geographical area represented by theDCM.
 22. An apparatus comprising: memory; and a processor coupled to thememory and configured to: generate for a Global Navigation SatelliteSystem (GNSS) a differential correction map (DCM) representing anon-planar surface of differential corrections that varies across ageographical area represented by the DCM, the differential correctionsbeing based on a reference station constellation of GNSS referencestations having respective locations spanning the geographical area, theGNSS reference stations including physical or virtual referencestations.
 23. The apparatus of claim 22, wherein the reference stationconstellation includes a virtual reference station constellation (VRSC)including virtual reference stations approximately uniformly spacedacross the geographical area represented by the DCM.
 24. The apparatusof claim 23, wherein the processor is further configured to: select theVRSC from a set of predetermined VRSCs, wherein the processor isconfigured to generate by generating the DCM for the selectedpredetermined VRSC.
 25. The apparatus of claim 24, wherein the DCMincludes: GNSS time; and a location of only one of the referencestations; an identifier of the VRSC; an identifier of each GNSSsatellite in view of the geographical area represented by the DCM; andthe differential correction for each reference station for eachfrequency of each GNSS satellite in view.
 26. The apparatus of claim 22,wherein the reference stations are non-uniformly spaced across thegeographical area represented by the DCM.
 27. The apparatus of claim 22,wherein the DCM includes: GNSS time; locations of the referencestations; an identifier of each GNSS satellite in view within thegeographical area represented by the DCM; and the differentialcorrection for each reference station for each frequency of each GNSSsatellite in view.
 28. The apparatus of claim 22, wherein the processoris configured to generate by generating the DCM from a two-dimensionalpolynomial surface fit to the differential correction associated witheach frequency and each GNSS satellite in view of the geographical arearepresented by the DCM.
 29. The apparatus of claim 28, wherein the DCMincludes: GNSS time; an origin of the two-dimensional polynomialsurface; an identifier of each of the GNSS satellites in view of thegeographical area represented by the DCM; and polynomial coefficientsassociated with the two-dimensional polynomial surface.
 30. Theapparatus of claim 22, wherein the processor is further configured toprovide the DCM to a transmitter configured to transmit the DCM over abroadcast coverage area via a radio signal.
 31. The apparatus of claim22, further comprising a network interface unit coupled to the memoryand the processor and configured to communicate with a communicationnetwork, such that the DCM is accessible over the communication network.32. The apparatus of claim 22, wherein the processor is furtherconfigured to: access from a database aggregated differential correctioninformation, wherein the processor is configured to generate bygenerating the DCM from the aggregated differential correctioninformation.
 33. The apparatus of claim 22, further comprising: atransmitter, wherein the processor is further configured to generatemultiple DCMs based on multiple reference station constellations, andthe transmitter is configured to alternately transmit the multiple DCMs.